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Abstract 

Equations describing the resistive state of a layered superconductor with 

anisotropic pairing are derived. The similarity with a stack of Josephson 

junctions is found at small voltages only, when current density in the direction 

perpendicular to the layers can be interpreted as a sum of the Josephson 

superconducting, the Ohmic dissipative and the interference currents. In the 

spatially uniform state differential conductivity at higher voltages becomes 

negative. Nonuniformity of the current distribution generates the branch 

imbalance and violates the Josephson relation between frequency and voltage. 
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Layered high-Tc superconductors are known to exhibit intrinsic Josephson effect when 
the current is ffowing across the conducting layers (see e.g. |I| and recent measurements 0). 
Though such materials can be considered as a stack of 2D-superconductors linked by the 
Josephson coupling 0,^, the theory of the ac Josephson effect in tunnel junctions cannot 
be directly applied to describe the resistive state of layered superconductors. In a system 
of series connected junctions the electric field is located mainly in the insulating barriers 
due to screening by the electrons in the metal, and the superconducting banks are in the 
equilibrium state. This results in many important consequences including the Josephson 
relation between voltage and frequency. But in the layered materials the superconducting 
layers are of atomic thicknesses and one must not ignore nonequilibrium effects which are 
related to perturbations of the quasiparticle distribution in the superconductor. On the 
other hand, many experimental evidences for (i-wave or nearly (i-wave symmetry of the 
superconducting order parameter in layered high-Tc superconductors were given last years, 
and a compatibility of the experimental data with d-wave scenario was shown in many 
theoretical works (see e.g. |]^,^] and references therein). In this case the superconducting 
order parameter has nodes, i.e. the quasiparticle density is never exponentially small and 
the nonequilibrium effects due to quasiparticles become especially important. Thus, to 
understand the intrinsic Josephson effect in high-Tc superconductors one must take into 
account the nonequilibrium distribution of the quasiparticles and the relaxation processes 
in the resistive state. 

In this study we calculate current and charge densities in superconductors with 
anisotropic pairing as functions of the phase differences of the order parameter in neigh- 
boring layers and of the nonequilibrium scalar potential related to the quasiparticle branch 
imbalance, i.e. the difference between densities of electron-like and hole-like quasiparti- 
cles 0,1^. We find that the direct analogy with a stack of Josephson junctions is limited 
by nonequilibrium effects and scattering processes, the difference with Josephson junctions 
being the most pronounced at lower temperatures. This results in the negative differential 
conductivity and in the violation of the Josephson relation. Similar effects are expected in 
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layered superconductors with isotropic pairing too, but in the latter case the quasiparticle 
density drops down exponentially when temperature decreases and generation of branch im- 
balance violating the Josephson relation would be negligible in the most interesting region 
of low temperatures. 

In our calculations we use the quasiclassical theory of nonequilibrium superconductivity 
IP modified for the case of layered superconductors [10|JTT[1 . We solve the equations in the 
discrete Wannier representation for the Keldysh |]12| matrix propagator, Gnm- Its diagonal 
components are the retarded and advanced Green's functions, and g^, and its upper 
off-diagonal component, g^ , is related to the electron distribution function. We consider 
the hopping conductivity regime between the layers, t^r -C h, which corresponds to the 
case of Josephson interlayer coupling. Here tj_ is the overlap integral describing the elec- 
tron spectrum in perpendicular direction, e± = 2tj_cos dk±, d is the lattice constant in the 
perpendicular direction, and r is the momentum scattering time along the layers. This ap- 
proach bears some similarity to the interlayer diffusion model [|T3| in which the interlayer 
coupling is mediated through incoherent hopping processes, being neglected. We assume 
that a symmetry of the superconducting order parameter is imposed by the symmetry of the 
coupling potential in the selfconsistency condition: thus we do not address the question of 
the microscopic nature of the interaction resulting in such a symmetry. Then the equation 
for Gnm has the form 



where is the angle of the in-plane electron momentum, (■ ■ ■) means averaging over 0, and 
hn = ~'^cry^n(0) + yUn + o"2vp„. Furthermore, /x„ = {h/2)dxn/dt + e^n is the gauge invariant 
scalar potential in layer n, Xn is the order parameter phase in the n-th layer, is the electric 
potential, pn = (/i/2)Vxn— (e/c)A„ is the superconducting momentum parallel to the layers. 
An is the vector potential, and Anm = cos [{xn — Xm)/2] + i<Jz sin — Xm)/'^]- The Pauli 
matrices introduced above act on the spin indices of Green's functions. Each Gnm depends 
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on two times (or on two energies in Fourier representation), and on 0. 

The r.h.s. of (0) is an elastic-collision integral in Born approximation. Using Born 
approximation we neglect low-energy quasiparticle bound states created by impurities (see 
and references therein), and, hence, our results are applicable provided typical energies 
of quasiparticles are larger than the bandwidth of the impurity induced bound states, T > 



m/r. 

We calculate current and charge densities assuming the clean limit, TcT ^ h, and neglect, 
where possible, the pairbreaking due to elastic scattering. We assume also the case of 
frequencies hu much smaller than typical electronic energies: temperature and the amplitude 
of the gap, Aq. 

Consider a current flowing between the layers of a quasi two-dimensional superconductor, 
t± ^ Aq. When the current exceeds its critical value it produces nonequilibrium pertur- 
bations the scale of which is determined by the interlayer coupling, and for small t±, the 
value of the nonequilibrium potential fi is small as well. Then we may solve equations (|1]) 
perturbatively, considering tj_ and /i as small values. We calculate the off-diagonal in layer 
number components of in the linear approximation in t±, these components are related 
to the current density in the direction perpendicular to the layers. The diagonal in layer 
numbers component of Trg^ determines perturbations of the charge density, we calculate it 
up to the second order in t_i_ neglecting n in comparison to A. The structure of the solution 
can be demonstrated by the combination gj = Ann-ign-in ~ 9nn-i^n-in the integral of 
which over energies and yields the current density across the layers [0. In the Fourier 
representation we get 



e / A2 AM ^ dnp ^it±9i\e\ - \A\) 



g, = -^t^tanh— — — - — — ]S^-J2 



^2(/in-l -/^n)^2 / ,c C )1 {2) 

where S, s, and c are Fourier components of simpn, siny9„/2, and cosv9„/2, respectively, 
V^n = Xn — Xn-i IS the phasc difference between the layers. The first term in (g) describes 
perturbations of the retarded and advanced propagators, it is related to the supercurrent. 



_ -|- ^g~)2 _ I Ap. The last terms describe perturbations of the distribution func- 

tion and contribute to the quasiparticle current, np is Fermi distribution function, and 



a = e/^, ^ = ye^ — |Ap. Furthermore, z/ = (l/r)[(a) — (A/e)(A/,^)] is effective momentum 
scattering rate for quasiparticles. Finally, Ub = —^^^ — {^/^){^/0] is the effective 

branch imbalance relaxation rate. The branch imbalance is related to the nonequilibrium 
potential /i i7|j8|JT5[| . In the case of isotropic pairing elastic scattering does not contribute to 



the relaxation of the branch imbalance, and the latter relaxes via energy scattering processes. 
We shall consider the opposite case of the order parameter close to ci-wave symmetry, when 
the gap has nodes: (A(0))^ ^ (A(0)^). In this limit we may simplify the equations because 
both z/ and Uh do not depend on (p, u = (l/r)(a) and Ub = (l/r)(A^(0)a/e^). Then the 
expression for the current density between layers n and n — 1 acquires the form 

1 ft dti r/ /■*! dt2 



j± = jc sm ipn+ , , 
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where generalized conductivities are given by 

ab{t,h,t2) = a^x(^e-'^(*-*^)e-'^^(*-*^)), (5) 
Averaging over the angles and quasiparticle energies in (H-BI) is performed according to 



(■■■), = -/ de{9{\e\-\Am^{---)). 



de 

The first term in (^) describes the Josephson current, jc = hc^ / {ATie\\d) being the critical 
current, is the penetration length for a superconducting current perpendicular to the 
layers. The quasiparticle contribution is given by the last terms containing retardation ef- 
fects related to the momentum and branch imbalance relaxation. These effects correspond 
to factors {uj + iv) in (|^). In the limit of A ^ these terms reduce to the Ohmic cur- 
rent aN±E/{l + iur), where a^i. is the static conductivity in the normal state. In the low 
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frequency limit the quasiparticle contribution may be interpreted as the Ohmic and inter- 
ference current (the last term), the Ohmic current consisting of two contributions, one of 
which is related to the branch imbalance and to its relaxation. When is a slowly varying 
function of time, u; -C z>, the retardation effects can be neglected, and the terms depending 
on the phase difference acquire a simple form typical for the Josephson tunnel junctions: 

= sm ^„ + — + a. ^ cos , (6) 



where conductivities are given by with the exponents integrated over time ti, which 
results in factors 1/z/. Near one gets a = (Tn^. At low temperatures, T <^ Aq, the 
quasiparticle conductivity is of the order of the normal state conductivity, because the 
decrease of the normal carrier density upon cooling is compensated by the decrease of the 
effective scattering rate of quasiparticles by the same factor oc T/Aq. For the simplest 
angular dependence of the gap parameter with the d-wave symmetry, A = Aq cos 20, we 
find cr = 3crAr^/4. For the interference term we get CTj = {/S./T) a^i. at T ^ A, and 
(Tj ~ 0"Ar^/4 at lower temperatures. 

The value of n in the current density (|^) must be determined from the separate equation 
describing the branch imbalance dynamics. Such an equation is given by the Poisson's 
equation with the charge density calculated from the integral of Tr g^^ over energies. In the 
Fourier representation we get for the charge density p„ in the n-th layer 

- iujpn = -iu;7— /i„ - V(o-2||V/i„/e + zcjai||P„) + — J^i'^n - Jn-i), (7) 

where J„ = J2uj2(^2±{n]UJ,uJi,uj2)[fJ'n{i^2) - yUn-i(t^2)] - (Ti±{n;uJi,uj)2ihuJiipn{uJi) describes 
the flow of quasiparticles between layers n and n—1 generating the branch imbalance. Note 
that (j{n) in the equations for charge density depend on the branch imbalance scattering 
rate, they are different /^from the conductivities for current densities. 

• 1 —'?h h 

7=1- 7 -— CTfell = CTtvII -7 --rrr-, , . (8j 

[ijj + ll'i,) t[ijJ + IVj^ljJ + lV}yY 
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a2±[n) = aN±{ -, — —^7 )e- (10) 

t{uJi + tU){uJ + tUb){uJ2 + IVb) 

^From equation (|^ one can see that the branch imbalance 7^ is generated by a nonuni- 
form quasiparticle current flow. 

In the limit of small phase differences between the layers equations and ([7[-pi]D 

describe the linear response characterized by different conductivities for the response to the 
solenoidal and to the potential electric fields |jT6[ . 



Using the definition of /i„ we get 

T/ ^^^^ fn\ 

eVn - --g^ = fin - fJ'n-l, (H) 

where = $„ — $„„i is the difference of electric potentials per a layer. Thus, the r.h.s. of 
(prp describes violation of the Josephson relation. 

We consider, first, the uniform case. In the case of spatially uniform current distribution 
both VPn = and J„ = Jn-i, and using in (0) p = we find pn = 0. Then, according to 
(pA|), the Josephson relation between the frequency and the electric potential difference per a 
layer, Vn, is satisfied. In the limit of small frequencies the current-phase relation has the 
form similar to Josephson tunnel junctions, and the effects typical for Josephson junctions 
must be observed. However, the analogy is limited by the region of voltages and frequencies 
smaller, than the effective momentum scattering rate of the quasiparticles, Uqp ^ z/(e = T). 
At higher voltages the finite scattering time effects, which were neglected in (||), destroy the 
analogy. 

In principle, for typical parameters of a superconductor both current biased and voltage 
biased regimes in the resistive state are possible. For simplicity we shall concentrate on the 
case of the voltage bias, which can be realized, for example, when a capacity is connected in 
parallel with the superconductor. Then we can easily find the time dependence of the phase 
difference using the Josephson relation, and calculate the current density /^from equation 



3± = (^N±{7T77Z7Z-rVT?I^iH^ + + 7T^V4a;2 + i>2 cos cut]),, (12) 



where huj = 2eV and we omitted the layer index in the uniform state. Equation ([l2D 
resembles expression for the current density across a tunnel junction. It contains dc and 
ac components, the characteristic frequency and voltage of such a junction, huc = 2eVc, 
are determined by the momentum scattering time: Vc ~ (Ao/T)^(/i/er) at T ^ Aq, and 
Vc ~ (T/Ao)(/i/er) at T ^ Aq. At eV ^ hugp the dc current decreases with voltage 
increasing. Note, that at lower temperatures eVc ~ fii^qp, thus, the I — V curve starts to 
decrease with voltage already at V ~ Vc- Though numerically the result depends on the 
details of the angle dependence of the order parameter, qualitatively it can be illustrated by 
the explicit calculations with A = Aq cos 20 at temperatures T ^ Aq. 



The origin of the negative differential conductivity at high voltages can be interpreted 
as the decrease of the dissipation at frequencies higher than the quasiparticle scattering 
rate, u ^ Ugp, or, equivalently, using the analogy to the negative differential conductivity in 



miniband and the momentum scattering rate. In the latter case the chemical potential in 
the adjacent layers is shifted by the value of the voltage exceeding the width of the band of 
the allowed electronic states, therefore, the electron's energy in one layer corresponds to the 
forbidden states in the neighboring layer. So, the effect must be present in the normal state 
of layered conductors as well (see also [|18| ). 

The negative differential conductivity indicates to an instability of the uniform resistive 
state at voltages eV > huqp. The instability must result in a nonuniform current distribution 
in which, according to (|^ /i 7^ and the Josephson relation is violated. Nonuniform current 
distribution may be created also at lower voltages due to many other reasons e.g. due 
to Meisner effect, due to contacts, or due to nonuniformities of the material. Thus, the 
correction to the Josephson relation in (^) depends on the experimental conditions. For 
illustration we estimate such a correction for an artificial but easily treatable model of 
the nonuniformity created by a layer dependent impurity scattering time, r„. We consider 




(13) 



semiconductor superlattices [|r^ at voltages per period higher, than both the width of the 



8 



the case when variations of with n are small, and the nonequilibrium potential /i can be 
calculated perturbatively. Then using condition p„ = we get for the case of low frequencies 



Estimating the integrals in (|T^ for T ^ Aq we get from ( ITll) [Vn — {h/2e)ipn]/Vn ~ 
h~'^t\Tn{Tn+i — r„_i)(Ao/T) lu Aqt. Ouc cau see that the significant violation of the Joseph- 
son relation may be created even by small variations of r in different layers. 

Thus, generally speaking, the Josephson relation in layered superconductors with Joseph- 
son interlayer coupling and anisotropic pairing is satisfied under the special conditions of 
the uniform current distribution only, which is difficult to satisfy. Even in the ideally uni- 
form samples the uniform state is expected to be unstable at higher voltages because of 
the negative differential conductivity. This may be the origin of the irregular character of 
the I — V curves observed usually in the measurements of the intrinsic Josephson effect in 
high-Tc materials. 
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